Bound states in open coupled asymmetrical waveguides and quantum wires 
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The behavior of bound states in asymmetric cross, T and L shaped configurations is considered. 
Because of the symmetries of the wavefunctions, the analysis can be reduced to the case of an 
electron localized at the intersection of two orthogonal crossed wires of different width. Numerical 
calculations show that the fundamental mode of this system remains bound for the widths that we 
have been able to study directly; moreover, the extrapolation of the results obtained for finite widths 
suggests that this state remains bound even when the width of one arm becomes infinitesimal. We 
provide a qualitative argument which explains this behavior and that can be generalized to the lowest 
energy states in each symmetry class. In the case of odd-odd states of the cross we find that the 
lowest mode is bounded when the width of the two arms is the same and stays bound up to a critical 
value of the ratio between the widths; in the case of the even-odd states we find that the lowest 
mode is unbound up to a critical value of the ratio between the widths. Our qualitative arguments 
suggest that the bound state survives as the width of the vertical arm becomes infinitesimal. 

PACS numbers: 03.65.Ge,73.21.Hb,73.20.-r 



I. INTRODUCTION 

Trapped waves in open geometries have been known 
for sixty years nowi. Ursell^, in a theoretical and ex- 
perimental study of "beach" waves in a particular semi- 
infinite canal, found that the system has discrete, con- 
tinuous and mixed spectra, and predicted the existence 
of a confined resonance at a discrete frequency, while 
at a cutoff frequency the resonance extends a long way 
down the canal. These waves can be described by the 
Helmholtz equation subject to Dirichlet boundary con- 
ditions and therefore similar confinement properties are 
expected in the solutions of analogous acoustic or elec- 
tromagnetic problems. Indeed, about fifteen years after 
Ursell, long-lived resonance modes in open laser systems 
were reported^. Even more, since this setup is equivalent 
to that of electrons in open configurations described by 
the time-independent Schrodinger equation, by 1984, sev- 
eral studies already suggested that confinement is also a 
feature of the two-dimensional transport of charge carri- 
ers in ultrafine metal and semiconductor devices^. From 
then on, many theoretical and experimental studies have 
shown that, generally, by bending and crossing waveg- 
uides and quantum wires, bound (confined, trapped) 
states can be obtained with energies below the contin- 
uous spectrum^. Moreover, many similar quantum wire 
systems have been found with bound states embedded 
into the continuum^. 

Bound states just below the continuum may strongly 
influence the scattering of waves and charge carriers. 
This allows to detect them and it is also the source of 
new phenomena and applications. For instance, meso- 
scopic systems with similar geometries are common in 
nanoelectronics and it has been shown experimentally 
that they exhibit suppression (quenching) of the Hall 
resistance^, and enhanced bend resistance^. On the other 
hand, if the above mentioned bound states exist only due 
to a particular symmetry, when this symmetry is broken. 



resonances could be expected to arise. Therefore, the 
confinement of electromagnetic radiation at the intersec- 
tion of coupled waveguides opens a window for a new 
kind of resonator, which combines a very simplified spec- 
trum to high performances^. Similar phenomena should 
be also relevant for the performance of novel kinds of 
(cold) aton*iS, phonoi^ and laser waveguidesi^. Even 
more, the possibility of geometrically inducing confined 
states allows for the design of unused ways to create 
a Bose-Einstein condensate in a quasi-one dimensional 
Bose ga&i^. Last but not least, besides the possible prac- 
tical applications, it should be noted that such systems 
constitute a very interesting theoretical framework for 
studying the correspondence between classical and quan- 
tum dynamics. Note that, in this kind of arrangements, 
classical motion of point particles with finite energy is 
typically unbound because there are not forbidden re- 
gions along the wires. Therefore, a phase-space semiclas- 
sical approximation to dynamics seems to be not valid 
here, and this is particularly relevant if the correspond- 
ing classical motion exhibits chao&i^. 

A relevant question, from both the theoretical and ex- 
perimental point of view, is how strongly depend the 
properties of the bound states on the perturbation of the 
structure and, this way, on the mechanical imperfections. 
Amongst other factors, is perhaps the symmetry of the 
configuration what plays the fundamental role since it 
controls the coupling with the propagating modes. Some 
studies suggested that the isolated bounded modes are 
not affected by the level of symmetry and thus, they can 
exist under very general conditions'^. Nevertheless, a 
comprehensive analysis of the role of symmetry is still 
lacking. 

In a classical paper written more than twenty years 
ago Schult, Ravenhall and Wyldi^ studied the problem 
of an electron which can freely move in a two dimensional 
symmetric region of orthogonally crossed wires of finite 
width and infinite length. They showed that the ground 
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state of this system corresponds to a state where the 
electron is confined in the region of the crossing, with an 
energy which fahs below the threshold of the continuum, 
i?TH- Using two different numerical approaches (finite 
differences and a mode expansion) they estimated the 
ratio between the ground state energy and the threshold 
energy, £'i/i?TH ~ 0.66. They also identified a second 
bound state, with odd-odd symmetry, corresponding to 
a ratio E2/Eh ~ 3.72. The energy of this second state 
falls below the threshold of the continuum for odd-odd 
states. These results were later confirmed by Avishai et 
al Ji. More recently, Amore, Fernandez and Rodrigueai^, 
have used the Conformal Collocation method (CCM)i^ 
to obtain a precise value for the ground state of the in- 
finite cross, using a non uniform grid. They obtained 
Ei/Etyi « 0.659611. Trefcthcn and Betcko^a, on the 
other hand have obtained a precise value for the second 
bound state, i?2/£^TH ~ 3.71648 (notice that these au- 
thors report the value of E2 , from which the ratio can be 
obtained) . 

Different geometries are closely related to the symmet- 
ric cross. Bulgakov and collaborators have studied the 
properties of a configuration of non-orthogonal ( "scissor- 
shaped") crossed wires, observing the emergence of mul- 
tiple bound states below the continuum as the angle be- 
tween the arms is reduced^i. A particular case of T- 
shaped waveguide is obtained from the symmetric cross 
by desymmetrizing the region for even-odd modes, al- 
though this particular system does not support bound 
states. Nazaro\4^, analyzed the bound states of a similar 
configuration, but with arms of different widths. 

Taking into account the above mentioned relevance of 
the symmetries for the spectrum, in this paper we study 
a different modification of the case discussed by Schult 
et alii^, considering an asymmetric cross, with orthog- 
onal arms of different width. Our goal is to investigate 
the behavior of the bound states of this system, as one 
of the arms is enlarged, and the possibility of the disap- 
pearance of these states or the appearance of new bound 
states. The related L and T shaped configurations arc 
also analyzed. 

The paper is organized as follows: in Section we de- 
scribe the system of the asymmetric cross and we make 
qualitative predictions on the behavior of the bound 
states as one of the arms is enlarged; in Section IIIII we 
present the numerical results, which have been obtained 
using a collocation method; finally, in Section lTVl we draw 
our conclusions. 



II. THE ASYMMETRIC CROSS 

As we have mentioned in the Introduction, Schult, 
Ravenhall and Wyldi^ have proved that an electron 
which moves freely in an infinite symmetric cross is local- 
ized in the central region of the cross when it finds itself 
in the ground state. They also discovered a second local- 
ized state, which is the lowest energy state with odd-odd 



symmetry. 

The goal of our paper is to investigate the behavior 
of these modes (as well as the possible appearance of 
further localized states), as the width of one of the two 
arms is changed, thus obtaining the asymmetric cross, 
shown in Fig. [TJ In particular, we ask ourselves whether 



FIG. 1. (color online) Asymmetric cross 



the localized modes of the symmetric cross can survive to 
arbitrary perturbations of this kind: what happens when 
one of the arms is much smaller than the other? Would 
we still have bound states? What does this imply for the 
L and T shaped configurations? 

To answer our questions we need to solve the scalar 
Helmholtz equation 
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A*„(x,?/) = En'Unix.y) , 



(1) 



where {x,y) G O = {|a;| < Wx/'2.,\y\ <Wy/2} and we 
are using h/m = 1. The wave functions obey Dirich- 
let boundary conditions on 917. We have called Wx and 
Wy the widths of the two arms, defining /3 = Wy/wx (for 
(3^1 one recovers the symmetric cross). 

Before trying to answer our question in a quantita- 
tive way, we may attack the problem qualitatively. Let 
us consider the case /3 ^ 1, meaning that the width of 
the horizontal arm is much larger than the width of the 
vertical arm, and assume that there is at least one lo- 
calized solution to eq.([T]). This situation is represented 
in figlH the grey regions in this plot represent the for- 
bidden region, where the eigenfunctions of eq. ([T]) must 
vanish identically; we assume that the dimensions of the 
red rectangle are large enough so that the localized wave 
functions are exponentially small outside the rectangle. 

Under these assumptions, one may attempt an approx- 
imate solution of the original problem inside the rect- 
angle, and thus obtain accurate approximations to the 
bound states of eq. ([T]), if they exist. 

The fundamental idea is to represent eq. as a 
Schrodinger equation 



-iA + V^(a:,y) 



■^n{x,y)^En-^n{x,y) , (2) 
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FIG. 2. (color online) Highly asymmetric cross: Wx <C Wy 



call (l): if the area of V is much larger than the area of 
57, a large number of basis functions will be needed in the 
calculation to achieve a suppression of the wavefunctions 
in the classically forbidden regions. 

Since we are interested in the limit Lx ^ Ly, we may 
use a simpler and more direct approach: provided that 
the rectangle V is very thin, the lowest modes of eq. ([5]) 
on V will be dominated by the longitudinal modes of the 
basis, with negligible contributions from the transverse 
excited modes. 

This observation justifies the possibility of using a por- 
tion of the Hilbcrt space of the problem, corresponding 
to the states y). 

In this limit the "high-energy" modes essentially de- 
couple from the problem and we may work with an ef- 
fective Hamiltonian, obtained by averaging the original 
Hamiltonian with respect to the mode Ciiu)- 
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V{x,y) 



, ix,y)£n 



(3) 



To allow numerical calculation the exact potential 
V{x,y) is replaced by an approximate potential V{x,y): 



V{x,y) 







{x,y) G Q 

{x,y) i n 



(4) 



which takes a large but finite value Vq in the forbidden 
region. If Vo is large enough (much larger than the en- 
ergies of the lowest states), the wave functions of these 
states arc exponentially suppressed in the classically for- 
bidden region (the grey region in the figure), and the 
contribution from this region becomes negligible. 

The advantage of reformulating the problem in this 
way should be clear, since we may now work on the entire 
rectangle, {x,y) G V = {\x\ < L^, \y\ < Ly}, and use the 
orthonormal basis 
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H,ff{x) = / dy (:,{y)H Ci{y) 
where the effective potential is 



K//(x)^ / dy e,{y)Vix,y) 

J —Lt, 



, l^^l < Wx/1 

'^VojtlnCiivfdy , \x\>Wx/2 



(5) 



With this qualitative argument we find that the ef- 
fective potential behaves as a square well of finite depth 
(which we will denote by A), whose width is precisely the 
width of the smaller arm, Wx (see the left plot of figlSj. 
It is well known that the ground state of this potential 
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FIG. 3. (color online) Elfective potentials corresponding to 
figs. [2] and E] (left) and to figsH and[6] (right) 



and nx,ny = 1,2,.... 

Kaufman, Kosztin and Schulteit^ have used this ap- 
proach to calculate approximations to the eigenvalues 
and eigenfunctions of the negative Laplacian on finite 
domains, applying the Rayleigh-Ritz method to the 
Schrodinger equation with V{x,y) — >■ V{x,y) on a rect- 
angular domain fully enclosing the region. 

The accuracy of this method clearly depends both on 
the finite value of the potential step Vq used in the nu- 
merical calculation and the ratio between the areas of V 
and the portion of area of 57 falling inside V (which we 



is always bound, independently of the width and depth 
of the well. On these grounds we expect that the funda- 
mental mode of the asymmetric cross of figEI will always 
be bound, regardless of the width of the two arms, as 
long as they are finite! 

We can also study the remaining lowest modes in each 
symmetry class using the same approach. In Figs|31[S]and 
[S]we display the desymmetrized regions corresponding to 
the symmetry classes of modes with odd-even, even-odd 
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and odd-odd respectively. Repeating the same steps fol- 
lowed for the full domain of figUl we may easily convince 
ourselves that the effective potential of Figs|4]and|6]cor- 5 W 

responds to a finite well potential with a infinite wall on 
one side (see the left plot of figE]), whereas the case of - 
Figl5]is analogous to the one of Figl2l — 
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FIG. 4. (color online) T-shaped configuration; desymmetrized 
region for odd-even modes of the asymmetric cross. 
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FIG. 5. (color online) T-shaped configuration; desymmetrized 
region for even-odd modes of the asymmetric cross. 
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FIG. 6. (color online) L-shaped configuration; desymmetrized 
region for odd-odd modes of the asymmetric cross. 



• the fundamental mode of a domain with the shape 
of an infinite L (which is the lowest energy odd- 
odd mode of the symmetric cross) is bounded, but 
becomes unbounded as the width of any of the two 
arms is increased, while keeping fixed the width of 
the other arm; 

• the fundamental mode of an infinite T shaped do- 
main (even-odd modes of the asymmetric cross), 
shown in FiglH is not bound for the symmetric 
cross, but it becomes bound as the width 0, 
but remains finite; 

• the fundamental mode of the "rotated" T of fig. |4] 
(odd-even modes of the asymmetric cross), is not 
bound either for the symmetric cross or as Wx ^ 0. 

Notice that the predictions regarding the T-shaped in- 
finite domain in FigIS] are in agreement with the findings 
of Nazaro\«2^. 

In the next Section we describe the numerical methods 
that we have used to test these predictions and we present 
the numerical results obtained with these methods. 



III. NUMERICAL RESULTS 



The simple quantum mechanical problem for the po- 
tential well of the right plot of figlS] docs not admit lo- 
calized solutions when the potential well is cither very 
narrow or very shallow. 

These simple arguments allow us to make the following 
predictions: 

• the fundamental mode of the symmetric cross stays 
always bounded, as the width of one of the two arms 
is increased arbitrarily, but remains finite; 



The correctness of the analysis in the previous section 
strongly depends on the assumption of the decoupling of 
the "high-energy" modes which leads to an effective one- 
dimensional problem. However, we are interested in any 
perturbation of the symmetry, not only the large ones 
given by /3 ^ 1. Thus, here we verify, by numerically 
solving eq. ([1]), that the results obtained in the previous 
section are correct, and, even more, that our predictions 
can be extrapolated to the whole range of possible values 
of /3. 
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We now illustrate the method that we have used in 
our calculations. Upon a simple rescaling of the y axis, 
defining y' — y//3 ix' ~ x), wc write eq.([T]) in the form 
of a modified Hclmholtz equation 
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ipn{x',y') = En^n{x' ,y') , (6) 



which allows us to work on an infinite symmetric cross 
of width with /3 playing a role of a parameter. 

As we are interested in the bound states of eq. ([6]), we 
may "cut" the cross at some large but finite distances 
Lx and Ly, by imposing Dirichlet boundary conditions 
at X = ±Lj^ and y = ^Ly. If and Ly are large enough 
the wave functions of the bound states are exponentially 
small at larger distances and the error introduced by this 
approximation is negligible compared to the error due 
to the discretization. In general, for /3 ^ 1, one may 
choose Lx 3> iy, as the probability density is mostly 
concentrated on the wider arm. The maximum size of 
Lx and Ly is essentially dictated by the total number of 
collocation points corresponding to this choice. 

On the other hand, excessively large values of Lx and 
Ly should be avoided, since they increase the number of 
collocation points and therefore the computational power 
needed in the calculation. 

To discretize eq. © we use a collocation approach 
based on "tent functions" (TF): on the interval \x\ < L 
we define the uniform grid, whose A'^ — 1 points are 
Xk = 2Lk/N^ where N is an even integer and k = 
-N/2 + 1, -N/2 + 2, . . . , N/2 - 1. The tent function 
peaked at the point Xk is defined as: 



4>k{x) 



Xk + 1- 





Xk-l < X < Xk 
Xk < X < Xk+1 
X < Xk-l,X > Xk+1 



(7) 



A function f{x) obeying Dirichlet boundary conditions 
at X = ±L {f{±L) = 0) can be interpolated using the 
TF as 

JV/2-1 
fc=-JV/2+l 

For a higher dimensional problem, the same procedure 
may be followed, using multidimensional TFs which are 
the direct product of the TFs along each orthogonal di- 
rection: 

^ki,k2,---Mai^l^^2,---,Xd) = (f>kiixi) (j)kAx2) ■■■ 4>kA^d) 

In principle the interpolation of a d-dimensional func- 
tion on a d-dimensional region \xi\ < Li, requires M = 

The eigenvalue equation ([S]) may then be converted to 
a matrix eigenvalue problem calculating the matrix ele- 
ments of the Hamiltonian operator between all the func- 
tions of a set: the size of the matrix depends on the 
total number of functions used in the calculation, which 



rapidly grows with the dimensionality of the problem. 
A drastic reduction of the number of functions however 
can be achieved by considering only functions which are 
peaked at points internal to the cross. 

A second observation concerns the choice of A'^^, i.e. 
the number of collocation points in each orthogonal di- 
rection: although in principle, Ni {i = 1, . . .d) can take 
any integer value, it is convenient to pick values of Ni 
which allow to sample exactly the border of the cross. 
The eigenvalues obtained for different grids sampling the 
border of the domain provide a monotonous sequence of 
values, which allows to obtain a good approximation via 
extrapolation^^. 

Before calculating the eigenvalues of the asymmetric 
cross, we have tested this collocation method on the 
symmetric cross, for which a precise result is available^^. 
The method used in that case is the conformal colloca- 
tion method (CCM)iS,. Working with grids with A^ = 



Nx 



Ny — 80, 120, . . . , 880 we have calculated the low- 



est eigenvalue of the corresponding collocation matrix, 
obtaining a monotonous sequence of values: the most ac- 
curate value, corresponding to the finest grid, provides a 
ratio Ei/E^YL = 0.66166, which is just 0.3% above the 
value reported by Amore et al We have also performed 
a least square fit of the monotonous sequence of values 
with the functional form oi + a2/N^ + as/N"^^ + a^/N^^, 
obtaining ai = 0.65955, which differs just for the 0.1% 
from the value previously reportecU^. This test makes us 
confident that our results are precise. 

We can now present the results obtained for the asym- 
metric cross, discussing separately the states with differ- 
ent symmetries. In all cases we expect that the wave 
function of a bound state will decay exponentially as one 
moves far away from the center of the cross; in particu- 
lar, the effective Hamiltonian approach that we have de- 
scribed in the previous section predicts \E'(x,y)|y gxod ~ 
^-x/e^ for |a;| ^ Wx/2, in the limit /3 1. We have 
found that there is at most only one bound state in each 
symmetry class. 



A. Even-even state 

We first study the lowest energy state with even-even 
symmetry of the asymmetric cross. In the third column 
of Tabic [T] we report the values of the ratio i?(°°)/i?TH, 
obtained for different values /?. We used three different 
sets, corresponding to L = 20 and A^ = 600 (set I), L — 
40 and A^ = 800 (set II) and L = 100 and N = 1600 (set 
III) . The fourth and fifth column report the values of the 
decay lengths £x and £y obtained fitting the numerical 
eigenfunctions with the asymptotic behaviors \I'(x, 0) ~ 
Q-x/£^ and '${0,y) w e~^/^», respectively for |x| ^ Wx/2 
and \y\ > Wy/2. 

These values are plotted in FigsI7]and[51 
The dependence on /3 is described very accurately by 
simple least square fits. In particular, within the accu- 
racy of our calculations, the fit of / Eth is consistent 
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with the survival of the bound state for /3 ^ oo. 

Notice that £x grows roughly cubically in /3 for /3 » 1, 
limiting the range of values of /? where the numerical cal- 
culation can be performed. As an example, in figslHland 
llOl we show different representations of our numerical re- 
sults for the wave function of the lowest even-even state 
for /3 = 2 (in these figures, as well in the other figures 
representing the wavefunctions, the region of the plot is 
smaller than the actual region where the numerical cal- 
culations were performed). 



B. Odd-odd state 

As it was found by Schult, Ravenhall and Wyldi^ the 
symmetric cross (/3 = 1) has a second bound state, which 
has odd-odd symmetry and therefore it is also an eigen- 
state of an infinite L. According to our qualitative discus- 
sion based on an effective Hamiltonian, we expect that 
this state may become unbound at some finite /?. 

In Table |lT] we report the results for the ratio 
^(oo) ^ jTjr^^ (second column) and for the longitudinal and 
transverse length scales, and £y (third and fourth 
columns). In this case the length scales are obtained fit- 



TABLE I. The ratio B^'^'V-Eth for the lowest even-even state 
of the asymmetric cross. Set I corresponds to using L = 20 
and TV — 600; Set II corresponds to using L = 40 and N — 
800; Set III corresponds to using L = 100 and iV = 1600. 





Set 


/-Bth 




£y 


1.0 


I 


0.662960 


1.098 


1.098 


1.1 


I 


0.723925 


1.335 


1.006 


1.2 


I 


0.774665 


1.613 


0.938 


1.3 


I 


0.816242 


1.936 


0.887 


1.4 


I 


0.849968 


2.309 


0.847 


1.5 


II 


0.879058 


2.770 


0.818 


1.6 


II 


0.900702 


3.267 


0.793 


1.7 


II 


0.918059 


3.830 


0.773 


1.8 


II 


0.931999 


4.463 


0.757 


1.9 


II 


0.943228 


5.172 


0.744 


2.0 


II 


0.952308 


5.960 


0.733 


2.1 


II 


0.959682 


6.832 


0.723 


2.2 


III 


0.965821 


7.959 


0.717 


2.3 


III 


0.970627 


9.053 


0.711 


2.4 


III 


0.974578 


10.252 


0.705 


2.5 


III 


0.977844 


11.564 


0.700 


2.6 


III 


0.980557 


12.992 


0.696 


2.7 


III 


0.982821 


14.542 


0.695 


2.8 


III 


0.984719 


16.219 


0.690 


2.9 


III 


0.986319 


18.026 


0.689 


3.0 


III 


0.987674 


19.965 


0.685 




FIG. 7. (color online) Ratio E'^^^'^/Eth as a function of /3 
for the asymmetric cross; the sohd hne is the least squares fit 
Ei/Eth\^^'^ = 1.003 - 0.505^-='-^^''2'^ + 0.165^-'^ ''''^ 




FIG. 8. (color onhne) £x and £y as a function of /3 for the 
lowest even-even state of the asymmetric cross; the solid 
and dashed hnes are the least squares fits £^f = 0.472 + 
0.644 /3^-^°^ and £^1™ = 0.656 + O.USp-^-^"^ . 



ting the asymptotic behavior of the wave function with 
^(a;,Wy/3) w e~^'^^^ and ^{wy/3,y) w e"*'/^", respec- 
tively for I a; I ^ Wx/^ and \y\ ^ Wy/2, since the wave 
function vanishes for x ^ or y ~ 0. 

In FiglTT]wc report the results for E'^""^ / Eth of the 
Table |TT] and compare them with the least squares fit 



e[°°^/e 



TH 



FIT 



-15.61 - 15.56/3 



1.974 



34, 



0.9869 



Notice that the fit has a maximum corresponding to (3 — 
1.12288. 

In Fig[T2] wc display the length scales £x and £y as 
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FIG. 9. (color online) Wave function of the lowest even-even 
state for (5 = 2. 
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FIG. 10. (color online) Wave function of the lowest even-even 
state for /3 = 2. 



functions of /3, and compare them with the fits 
,FiT 0.0108614 



L 



1 - 0.992232/30 0673327 



and 



I FIT 



0.514767^ 



0.81312 

«12.3632 • 



We believe that this plot offers convincing evidence of the 
existence of a sing ularity in 4 close to /3i°°^ = 1.2279: 



TABLE II. The ratio E^°°'> /£th for the lowest odd-odd state 
of the asymmetric cross. All the results are obtained using 
L = 100 and iV = 1600. 



p 


£;(°°V-Eth 




y 


1.00 


3.72042 


1.332 


1.332 


1.01 


3.75611 


1.465 


1.232 


1.02 


3.78877 


1.623 


1.148 


1.03 


3.81839 


1.816 


1.076 


1.04 


3.84499 


2.055 


1.014 


1.05 


3.86855 


2.359 


0.960 


1.06 


3.88909 


2.760 


0.912 


1.07 


3.90659 


3.313 


0.870 


1.08 


3.92107 


4.125 


0.832 


1.09 


3.93252 


5.429 


0.797 


1.10 


3.94095 


7.875 


0.766 


1.11 


3.94635 


14.119 


0.739 


1.111 


3.94673 


15.322 


0.735 


1.112 


3.94707 


16.745 


0.733 


1.113 


3.94739 


18.456 


0.730 


1.114 


3.94767 


20.549 


0.728 


1.115 


3.94792 


23.167 


0.726 


1.116 


3.94815 


26.528 


0.722 



4.00 
3.95 
3.90 
3.85 
^^3.80 
3.75 
3.70 
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FIG. 11. (color online) Ratio E'-°°'>/Eth as a function of /3 
for the asymmetric cross; the solid line is the least squares fit 

I FIT 

e[°°^/Eth\ = -15.61 - 15.56/3i'9^^ +34.88/30-9869. 



for /3 > the wave function becomes unbound, in 

agreement with our earlier prediction. As an example, 
in figsll3land ll4l wc show different representations of our 
numerical results for the wave function of the lowest odd- 
odd state for /3 = 1.1. 
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FIG. 12. (color online) and iy as a function of /3 
for the lowest odd-odd state of the asymmetric cross; the 
solid and dashed lines are the least squares fits ^j,!^'""" = 
0.0108614/(1 - 0.992232^''°''''^^2'') and lyf^'^ = 0.514767 + 
0.81312/3^^^-^®^^. The vertical line corresponds to the critical 
value = 1.2279 where l^^l^^"^ is singular. 




FIG. 13. (color online) Wave function of the lowest odd-odd 
state for P = 1.1. 




X 

FIG. 14. (color online) Wave function of the lowest odd-odd 
state for /3 = 1.1. 

unbound; a second region, where the wave function is 
bound, but mainly localized in the vertical arm (see for 
instances figsfT5land [T6| . and a third region, where the 




FIG. 15. (color online) Wave function of the lowest even-odd 
state for /3 = 1.55. 



C. Even-odd state 

In Table|III]we report the results for E^''°'> / E^b.-, and 
iy of this state, calculated using L = 100 and N = 1600. 
The behavior of the lowest even-odd state is rich and it is 
characterized by three different behaviors: a first region, 
below a critical value of B. where the wave function is 



wave function is still bound, but mainly localized in the 
horizontal arm (sec for instances figs ITTI and [T8| . 

The ratio i?(™'/i?TH is plotted in FiglTOland compared 
with the fit 



TH 



FIT 



3.96521 



11.0968 10.156 



/3' 



7.38923 



3.69462 ■ 
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TABLE III. The ratio E'^""^ /Eth for the lowest odd-odd state 
of the asymmetric cross. All the results are obtained using 
L = 100 and = 1600. 
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FIG. 16. (color online) Wave function of the lowest even-odd 
state for (3 = 1.55. 




FIG. 17. (color online) Wave function of the lowest even-odd 
state for /3 = 3. 



As for the even-even state, this behavior is consistent 
with the survival of the bound state for /3 — > oo. 

In Fig. [20] we plot and iy and compare with them 
the fits 

41™ 0.129087/3^-2^^^2 -t- 0.258401 

and 

FIT _ 0.604788 

*-v\ "j 2.16467 • 

The singularity of iyf^"^ is located at /3i''°'' = 1.513 (the 
vertical line in the plot) and represents the critical value 



where the bound state appears. Notice also that 4 1 ~ 
lyf^'^ for/3 = 2. 

D. Odd-even state 

In the case of the odd-even states we have not been 
able to find a bound state for any of the values of /3 
considered. Once again this behavior is consistent with 
our earlier predictions. 
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FIG. 18. (color online) Wave function of the lowest even-odd 
state for /3 = 3. 






FIG. 20. (color online) and iy as a function of /3 
for the lowest even-odd state of the asymmetric cross; 
the solid and dashed lines are the least squares fits 
lx\™ = 0.129087,g^-^^^" + 0.258401 and ^^1^'"^ = 
0.604788/(1 - 2.16467/3"^®''^*''). The vertical line corre- 
sponds to the critical value /3i°°' — 1.513 where ^yl'^'^ is 
singular. 
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FIG. 19. (color online) Ratio E'~''°' /Sth as a function of /3 
for the asymmetric cross; the solid line is the least squares fit 

, ^ I FIT 

7Ti(oo) ij-, _ o oKt^OI J_ 11.0968 10.156 

-til /-ciTH — ^.96521 -f g i Aim'i-i ~ g-i-mm-i ■ 



IV. CONCLUSIONS 

For years now, it has been verified theoretically, as well 
as experimentally, that the spectra of waves or quantum 
particles moving in several open two-dimensional systems 
could exhibit bound states. Since these bound states 
arise from the symmetry in the system, for applications 
it is important to find out how the perturbation of the 
given symmetry affects the corresponding spectrum. 

In this paper we have shown that the question about 
the role of the symmetry of several setups can be an- 
swered by studying the simple problem of an effective 



one-dimcnsional quantum well. We support this conclu- 
sion by means of very accurate numerical solutions of the 
full Helmholtz equation with Dirichlet boundary condi- 
tions. It is worth to notice that the conclusions drawn 
from the simple problem are valid even for values of the 
ratio of the width of the crossing branches (/3) for which 
the one-dimensional approximation not longer holds. 

In particular, we have shown that the bound state of 
the symmetric cross persists for any finite value of /3. A 
similar, but richer, behavior is obtained for the T-shaped 
configuration as the width of the vertical bar of the T 
tends to zero (i.e., — > oo): from ~ \ vrp to 0i, ~ 
1.513 the wave function is unbound; above that critical 
value the wave function becomes bounded, but mainly 
localized in the vertical bar of the T; finally, if is still 
increased, the wave function remains bound, but mainly 
localized in the horizontal bar of the T. Of course, in 
any of these two setups, the bound state must dissapear 
as soon as the width of the decreasing branch reaches 
zero. Therefore, this seems to imply that adding just a 
very small perturbation on a very large one-dimensional 
quantum wire or waveguide may have a striking impact 
on the transport of charge carriers or waves along it. 

The situation is different for the other two cases we 
analyzed. For the L-shaped configuration we have found 
that the bound state becomes unbounded at the critical 
value 0i, = 1.2279 of the ratio of the width of the inter- 
secting branches. Finally, for the T-shaped configuration 
no bound state arises when the width of the horizontal 
arm is decreased towards zero. 

How strong are these effects as to be observable and 
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appliable will depend on how strong is the departure of 
the actual experimental setup from the ideal configura- 
tions considered. Nevertheless, if the dependence on the 
symmetry found here is confirmed by future experiments, 
then it could be used as a switch-like mechanism for trap- 
ping and untrapping waves and charge carriers. This 



would give rise to a number of important applications. 
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